Response to targeted perturbations for random walks on networks 
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We introduce a general framework, applicable to a broad class of random walks on networks, that 
quantifies the response of the mean first-passage time to a target node to a local perturbation of 
the network, both in the context of attacks (damaged link) or strategies of transport enhancement 
(added link). This approach enables to determine explicitly the dependence of this response on 
geometric parameters (such as the network size and the localization of the perturbation) and on the 
intensity of the perturbation. In particular, it is showed that the relative variation of the MFPT 
is independent of the network size, and remains significant in the large size limit. Furthermore, in 
the case of non compact exploration of the network, it is found that a targeted perturbation keeps 
a substantial impact on transport properties for any localization of the damaged link. 

PACS numbers: 



Since the pioneering works of Erdos and Renyi on ran- 
dom graphs PP, complex networks have been extensively 
studied, either for their theoretical interest, or for their 
applications in various areas such as social and economi- 
cal sciences, biology or epidemiology (see for instance [2J- 
13] and references therein). Among other features, trans- 
port properties on networks play a crucial role, and in 
this context, random walks on networks have appeared 
as a prototypical example of dynamical process which has 
been extensively studied over the past few years [SHE] . In 
particular, the target search problem has raised much at- 
tention because of the variety of its applications [9] , and 
the mean first-passage time (MFPT) to a target node 
is the most commonly used quantitative indicator of the 
efficiency of a search process on a network [TUHH]- A 
general framework to calculate MFPTs on networks has 
been presented in [5] and since then explicit results have 
been obtained both for fixed [TT] and averaged positions 
[I2WT5] of the starting point. 

In parallel, increasing interest in modeling failures or 
attacks in networks has developed [T2], for instance in 
the context of epidemic spreading |17j . virus attack on 
networks [TH], or electrical blackout [TH1 [20]. Most of 
these studies deal with static properties of networks af- 
ter a given perturbation (typically the removal of nodes 
or links), with the notable exception of |21j . In particu- 
lar it has been showed that scale free networks are very 
resilient to random perturbations, while the targeted re- 
moval of a hub can have dramatic consequences. Very 
recently, it has also been put forward that in the case of 
interdependent networks a broader degree distribution 
increases the vulnerability to random failure [T5] . 

In this context, quantifying the impact of targeted per- 
turbations of a network on the search efficiency and more 
generally on transport properties is an important and yet 
unexplored problem. In this letter, we provide a general 
framework that quantifies the response of the MFPT to a 
target node to a local perturbation of the network, both 
in the context of attacks (damaged link) or strategies 
of transport enhancement (added link). This approach 
enables to determine explicitly the dependence of this re- 



sponse on geometric parameters (such as the network size 
and the position of the perturbation) and on the intensity 
of the perturbation. In particular, it reveals that the rela- 
tive variation of the MFPT is independent of the network 
size N in the large N limit, and remains significant even 
for very large networks ; additionally, in the case of non 
compact exploration of the network, a targeted perturba- 
tion keeps a substantial impact on transport properties 
for any localization of the damaged link. 

We consider a discrete time random walker on a net- 
work of N nodes, and assume that the transition proba- 
bilities R wu from u to w of the walker are such that an 
stationary distribution P x exists. In addition, we denote 
P™ y the probability to be at site x after n steps for a ran- 
dom walk starting from site y. In what follows, we will 
be interested in the influence of targeted perturbations 
of the network on the MFPT (Tts) from a starting site 
5 to a target site T. 

We first introduce the general method to study the in- 
fluence of a single link perturbation on the MFPT, and 
take the example of the weakening or removal of a link. 
More precisely, we assume that the transition probability 
R wu from u to w is changed by SR WU < 0. Without loss of 
generality, we here assume that this perturbation is com- 
pensated on the probability R uu to stay at u during the 
elementary step, ie 5R UU = —SR WU , and that all other 
transition rates are unchanged. Note that the important 
particular case of a broken link is then given by taking 
5 R W u — —Rwu- In the sequel, all quantities denoted with 
a prime correspond to the perturbed situation. 

The MFPT in the perturbed situation can be calcu- 
lated by first noting that the perturbation affects only 
the trajectories that pass through u before reaching T, 
so that, for any starting point x : 

Ft*) ~ ( T Tx) = PL ((Tk) ~ (Tm)) , (1) 

where P£ x is the splitting probability to reach u before 
T starting from x. Writing next the backward equation 
for the MFPT [22] : " 

{T^) = 1 + Y j Ku{Ttu) (2) 
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both in the perturbed and unperturbed situations, we 
obtain 

(T Tu ) - (Ttu) = « T tJ - {T Tv )) 

V 

+SR wu {(T Tw ) - (Ik». (3) 

Using Eqs ([!]), ([3]), the variation of the MFPT starting 
from site u is found to be given by : 

irpl \ irp v _ SR wu ((T Tw ) - (T Tu )) 

KJ-Tu) ~ \J-Tu) - _ pT SJf — n _ pT y W 

± J tm T "-"^uitiA 1 1 uw) 

where we have introduced PJ U = ^ v Rv U Puvi defined 
as the probability to come back to u before reaching T. 
Using (IT]), we finally obtain the relative MFPT variation 
S TS = (fe) - (T TS ))/(Tts) for any starting site S: 



Sts — 



P^s SR wu {(T Tw ) - (T Tu )) 
(T TS )l-pT u + 5R wu {l-PT w 



(5) 



Using [TTJ [53], Sts can be expressed as a function of 
the perturbation SR WU , the unperturbed stationary dis- 
tribution P x , and the pseudo-Green function H xy = 

oo 

y~^(P"^ — P x ) of the unperturbed problem, giving : 

n=0 
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Htt — Hts 
5R wu [Pt{H u s — H u t) 



P u (Htt — Hts)} 



Pt + SR W u[Pt(Huu — Huw) + Pu{Htw — H 



Tu) 



(6) 



This central result has several important implications. 

First, we stress that in the particular case of a regular 
d-dimensional hypercubic parallelepipedic network with 
constant probability transitions between nearest neigh- 
bors the pseudo Green functions are known exactly [2"3] . 
Using next that the stationary probability is in this case 
uniform (P x — 1/N for any node x), Eq.([6| provides an 
exact and fully explicit result for the effect of an arbitrary 
modification of a given link on the MFPT. 

Second, in the case of more general networks, but pos- 
sessing scale-invariant properties, the dependence on the 
geometrical parameters can still be determined, by tak- 
ing the large network size limit. More precisely, this limit 
can be conveniently discussed when the random walk is a 
scale-invariant process, ie when the infinite volume prop- 
agator satisfies P xy oc n~ dl / d ™Il{\x — y\/n 1 / dvj ) where d w 
and df denote respectively the walk dimension and the 
fractal dimension of the network, and where \x — y\ de- 
note the distance between nodes x and y. Indeed, in this 
case all the differences H xy — H xz that enters Eq. ^ 
can be rewritten in terms of differences H xx — H xy which 
satisfy in the large volume limit [TT] : 
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FIG. 1: Relative variation of the MFPT as a function of the 
distance \S — T\ between source and target (here the relative 
positions of T, u and w are fixed and averages over triplets 
are taken). Numerical simulations (symbols) are compared to 
the approximation (6[7l (plain lines) for 3 different network 
sizes. Circles stand for 3D regular lattice (d w — df = — 1), 
triangles for 2D critical percolation cluster (d w — df ~ 0.98), 
crosses for a (2,2) flower (d w — df). 



where A, B are numerical constants depending only on 
the infinite volume propagator. Equations ([6])-([7| have 
two consequences, (i) The independence of N of Eq.([7]) 
readily gives that the relative variation Sts is indepen- 
dent of N in the large volume limit. This quite unex- 
pected effect, illustrated by the data collapse for differ- 
ent volumes in Fig. [T] on various examples of networks, 
implies that the effect of a targeted perturbation is not 
diluted but remains finite even for extremely large net- 
works. Actually this universal asymptotic independence 
on N of the relative variation 5ts for fixed starting node 
S strongly differs from the relative variation 5t of the 
MFPT averaged over S defined by (T T ) = £ s P s (T T s)- 
More precisely, assuming next detailed balance, it can be 
shown using symmetry relations of H xy that 



(T T ) - (Tt) 
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(8) 



Pt + SR wu [Pt(H uu — H uw ) + P u (Ht w ~ Ht u ) 
In the large volume limit, this exact expression leads to 
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where we have derived the asymptotic expression of Htt 
using [T5l [25] . In other words, as for the relative vari- 
A + B\x — y\ d ™~ d f if d w < df ation St, the iV independence is recovered only in the 

A + B In \x — y\ if d w = df (7) case df > d w of so-called non compact exploration, while 
B\x — y\ d ™~ d i if d w > df a strong dependence on N is found in the opposite case of 
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compact exploration (d w > df). This effect is illustrated 
in Fig. [2j (ii) The asymptotic form Q used in Eq. ^ also 
provides the explicit dependence of Sts on the relative 
distances between the nodes S, T, u, w. Such dependence 
has been checked numerically (see Fig. [T]) for various 
networks such as regular euclidian lattices (d w =2), 2D 
critical percolation clusters (d w = 2.88 and dt = 91/48), 
and (u, u)-flowers which are recursive fractals defined in 
[2"rJ] . Fig. [I] reveals a very different behavior in the case 
of compact exploration (illustrated by critical percola- 
tion clusters) for which Sts vanishes at larges distances, 
and in the case of non-compact exploration (illustrated 
by 3D regular lattices), for which Sts remains finite even 
for very large distances. Noteworthily this shows that 
in the non compact case a targeted perturbation keeps a 
substantial impact on transport properties for any local- 
ization of the damaged link. 
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FIG. 2: Relative variation of the MFPT averaged over S, St 
for different network sizes N for ID (black), 2D (red) and 
3D (blue) cubic lattices and a 2D critical percolation cluster 
(green). The relative positions of T, u and w are fixed for all 
networks of a given kind. The circles stand for the simulated 
St, the black lines for the theoretical prediction where C 
is a fitting parameter. 

Third, in the general case of a random walk on an 
arbitrary network, where the pseudo-Green functions can 
be difficult to evaluate, Eq. Q still gives explicitly the 
functional dependence of Sts on the perturbation SR WU , 
which takes the form 

<5 TS _ DSR w / 10 \ 
E + 5R WU 

where D and E do not depend on SR WU (note that E does 
not depend on S either). This general form has been val- 
idated by numerical simulations in Fig. [3] Additionally, 
provided that the differences H xx — H xy involved in Eq. 
^ have a finite limit in the large volume regime, D and 
E turn out to be independent of N. This shows that the 
independence of Sts on A still holds in this case, and 
makes this property very robust. 
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FIG. 3: Relative variation of the MFPT as a function of the 
u — > w link perturbation, for given T, u and w sites, and two 
different S for each network : 2D critical percolation clusters 
(circles) and (3, 3)-flower of generation 3 (triangles). For the 
flower network, the perturbed link leads to T. Numerical 
simulations are fitted with Eq. 
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Finally, it should be noted that the relative variation 
Sts remains rather weak for an arbitrary perturbed link. 
We however stress that such local attack of a network is 
not affected by dilution effects and remains finite even in 
the large volume limit; additionally, in the non compact 
case, it is also widely independent of the localization of 
the perturbation and non vanishing even for a very re- 
mote damaged link. Furthermore, the effect of a local 
perturbation can become much stronger if targeted to a 
link directly leading to the target, as can be expected 
intuitively (see Fig. [3]). 

Importantly, the above formalism can be extended to 
tackle the reciprocal problem of enhancing instead of 
damaging the transport abilities of a network. As a first 
step in this direction, we quantify the effect of adding a 
new link between two nodes. The definition of the per- 
turbation has to be slightly modified in the case of an 
added link, since the case studied above is ill defined for 
SRwu > and R uu — 0. We assume now that there is 
initially no link between u and w {R W u = 0), and con- 
sider a perturbation SR WU > 0. In turn, we set for all 
neighbors vofu that SR VU = —SR wu /k(u), where k(u) is 
the initial connectivity of u (note that SR WU is assumed 
to be small enough so that all transition probabilities are 
positive). In this case, the equivalent of Eq. ([5| can be 
obtained along the same line and reads: 



Sts — 



P T 

r uS 



TS/ 



SR W u (\Ttw) 
1 - PL + SR W 



(T T u) + 1) 
( pT _ pT 



(11) 



which, as previously, can be expressed only in terms of 
SR WU and pseudo-Green functions. This theoretical pre- 
diction is plotted against numerical simulations on the 
example of parallelepipedic networks in Fig. [4] (black 
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FIG. 4: Relative variation of the MFPT in response to the 
addition of a new link, in the directed (red and black lines 
and symbols) and bidirectional case (green line and symbols) 
for 3D regular lattices. Simulations (symbols) are plotted 
against the theoretical prediction (Eq. (Ill and its extension, 
plain lines). Black circles: S, T, u and w are distinct, and 
the new link starts from a target neighbor. Red circles: the 
new link is between 5* and T. Green triangles : addition of a 
link between two 3D regular lattices (with 5R WU = SR um ), S 
being in the first lattice, and T in the second. Here we define 
Sst = (MFPT - min(MFPT)) / min(MFPT) . 



and red lines). Note that the MFPT can be decreased 
very significantly if the new link points to the target (red 
line of Fig. Last, we stress that the case of a single 
added bidirectional link connecting two nodes u and w of 
initially distinct networks M\ , A/2 can be obtained using 
the same method, and yields an explicit results for the 
MFPT which is displayed in Fig. [4] (green line). This 
constitutes a first step in designing interdependent net- 
works as introduced in |19j . 



To conclude, we have presented a general framework 
that quantifies the response of the MFPT to a target 
node to a local perturbation of the network, both in the 
context of attacks (damaged link) or strategies of trans- 
port enhancement (added link). This approach enables 
to determine explicitly the dependence of this response 
on geometric parameters (such as the network size and 
the position of the perturbation) and on the intensity of 
the perturbation. It reveals that the relative variation 
of the MFPT is independent of the network size, and 
remains significant even in the large size limit. Addition- 
ally, in the non compact case a targeted perturbation 
keeps a substantial impact on transport properties for 
any localization of the damaged link. 
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